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ABSTRACT 


An analysis of historical trends in diversification can suffer from a variety of defects and limitations. Among 
these is the potential failure of genus and family compilations to yield sufficient information on biological events, 
such as mass species extinctions. Another is the inadequacy of taxonomically based estimates of diversification 
to reveal patterns in changing community structure and morphology, which are important in evaluating adaptive 
evolution. These deficiencies are examined for the plant fossil record and by means of computer simulations of 
mass extinction events. Comparisons among species, genus, and family diversification patterns reveal qualitative 
differences sufficient to warrant caution in the use of taxonomic ranks above the species level. However, genus 
tabulations appear to be an adequate surrogate database for inferring patterns of species origination. Computer 
simulations revealing patterns of family extinctions may be inadequate for evaluating mass species extinction 
events since peaks in species and family extinctions rarely, if ever, coincide. A pluralistic approach to evaluating 
diversification is advocated involving the examination of biotic changes within assemblages and trends in 
morphological, anatomical, and reproductive evolution. 


The notion of diversity in paleontology 
differs little from its use in ecological studies 
(Raup & Stanley, 1971; Whittaker, 1977). 
Diversity can be defined as the number of 
taxa in a community or as the synthetic char- 
acteristic of taxonomic richness and equita- 
bility, i.e., the relative evenness of the im- 
portance values of taxa within a sample (Lloyd 
& Ghelardi, 1964). However, the nature of 
the fossil record precludes a direct compari- 
son of diversity between past and present 
biotas (Krassilov, 1975; Padian & Clemens, 
1985). In ecology, taxonomic richness re- 
flects the number of species in a sample of 
standard size, while measures of equitability 
involve some gauge of the abundance or pro- 
ductivity of one species divided by the im- 
portance values of all other species within the 
sample (Pielou, 1977). By contrast, fossil as- 
semblages usually reflect time-averaged sam- 
ples of a geographically ill-defined area. Fur- 
ther, terrestrial organisms are preserved where 
they were buried, not generally where they 
lived, and most often after significant periods 


of decomposition (cf. Niklas et al., 1980; Pa- 
dian & Clemens, 1985). Frequently, trans- 
port to sites of burial results in the disarticula- 
tion of land plants and animals, and in 
deposition of parts in sedimentologically dif- 
ferent microenvironments (cf. Kidwell, 1986). 
Consequently, there is great difficulty (1) in 
recognizing a species from its parts, (2) reas- 
sembling the components of a community, 
and thereby (3) arriving at a measure of di- 
versity that is comparable to contemporary 
ecological studies (Raup, 1976, 1979; Niklas 
et al., 1980; Knoll, 1984; Benton, 1985). 
Despite its limitations, the fossil record pro- 
vides a potentially valuable perspective on a 
number of evolutionary issues. The paleon- 
tologist can compare patterns of diversifica- 
tion over billions of years of Earth’s history 
and can track the origin, radiation, and even- 
tual taxonomic diminution of organisms that 
no longer exist or are rare in present biotas. 
Although the nature of paleontological data 
defines the temporal and taxonomic resolution 
with which evolutionary issues can be ad- 
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dressed, recognition of large-scale patterns 
that develop, change, and disappear over long 
intervals of time confers a perspective that 
augments neontological studies. 

This paper is concerned with changes in 
tracheophyte diversity over the last 400 mil- 
lion years of the Phanerozoic. The fossil re- 
cord of terrestrial vascular plants is reviewed 
and discussed within the context of large-scale 
patterns in composition and taxonomic rich- 
ness of assemblages through time. The bulk 
of the data to be discussed has been presented 
elsewhere in a variety of forms (Niklas et al., 
1980, 1985; Tiffney, 1981; Knoll, 1984, 
1986; Knoll et al., 1984). Consequently, 
many facets of tracheophyte diversification 
will be referenced but not discussed in detail. 
The principal focus of this paper will be on 
the limitations of the data and on the quali- 
tative statements about plant evolution that 
can be made from an inherently imperfect 
fossil record. 


LARGE-SCALE CHANGES AND 
TAXONOMIC LEVEL 


The species is the taxonomic level of choice 
for determining changes in diversity through 
geologic time. However, species-level com- 
pilations present a variety of procedural dif- 
ficulties, the most significant of which involve 
sampling errors (Signor & Lipps, 1982; see 
also Pease, 1985). Comparisons of global pat- 
terns in species diversification with higher 
taxonomic levels are desirable, because they 
provide insights into the information to be 
gained or lost by each type of compilation. 

The upper panel of Figure | shows a plot 
of the total species-richness of vascular plants 
through the Phanerozoic. The data are taken 
from primary paleobotanical citations tabu- 
lated by Niklas et al. (1985) and are segre- 
gated into three major categories of tracheo- 
phytes: pteridophytes, gymnosperms, and 
angiosperms. The plot shows a Silurian to mid- 
Devonian radiation of primitive pteridophytes 
followed by a Carboniferous proliferation of 
more advanced or derived pteridophytes and 
early gymnosperms. A Permo-Triassic de- 
crease in overall diversity precedes a more 
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or less constant level of species number until 
the middle of the Cretaceous when angio- 
sperms begin to radiate. Diversity of pteri- 
dophytes decreased during the Permo-Trias- 
sic and reached a relatively constant level 
throughout much of the Mesozoic and Ce- 
nozoic. With the advent and proliferation of 
angiosperms, the number of gymnosperm 
species decreased significantly through much 
of the late Cenozoic. 

The lower panel in Figure | shows large- 
scale temporal patterns in diversity at the 
family level. These data (from Knoll, 1984) 
were compiled independently from the data 
on species from Niklas et al. (1980). Data 
for pteridophytes and gymnosperms are not 
segregated, thus providing for direct com- 
parisons between species- and family-level di- 
versity for nonangiosperms and for angio- 
sperms. Differences in the diversity between 
Devonian and Carboniferous families are less 
pronounced than those seen in the top panel. 
As Knoll (1984) noted, the number of species 
in late Paleozoic families is higher than in 
early Paleozoic families, presumably due to 
the appearance of more complex or numerous 
morphological and anatomical features upon 
which species can be based. Ordinal changes 
in plant diversity compiled by Knoll (1984) 
conform in large measure to the trends seen 
in the family diversity plot. 

Despite the differences between species and 
family diversity, it is evident that the number 
of vascular plant taxa increased throughout 
much of the Paleozoic and has accelerated 
during the Cenozoic with the advent of an- 
glospermy. 

The value of comparing diversity among 
different taxonomic levels can be further il- 
lustrated by examining data for a specific 
group rather than the entire domain of tra- 
cheophytes. The selected group ought to have 
a fossil history that spans the Phanerozoic 
and thus obviates any idiosyncratic feature 
of a geologic period. Three major groups of 
vascular plants have a virtually continuous 
fossil record since the Upper Devonian: ly- 
copods, sphenopsids, and ferns (Stewart, 
1983). Of these three groups, the ferns are 
the most abundant through time. Hence, they 
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Large-scale changes in diversity of tracheophytes at the species- and family-level (upper and 
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lower panels, respectively) recorded per epoch through the Phanerozoic. Species global diversity data were 
taken from Niklas et al. (1980, 1985); family diversity data are redrawn from Knoll (1984). The time-scale is 
that of Harland et al. (1982): For the Carboniferous (C), T = Tournaisian, V = Visean, N = Namurian, W = 
Westphalian, S = Stephanian; for the Cretaceous (K), L = Berriasian-Barremian, A-A = Aptian-Albian, U = 


Cenomanian-Campanian, and M = Maastrichtian. 


are convenient for comparing differences in 
diversity based on compilations of species, 
genera, and families. The most complete 
treatment of fossil ferns is that of Boureau 
(1970). Unfortunately, this does not reflect 
recent advances in taxonomy nor does it pre- 
sent a critical evaluation of stratigraphic oc- 


currences (see for example Phillips, 1974; 
Stewart, 1983). Nonetheless, the data con- 
tained within this volume are useful because 
they are easily accessible and reflect the state 
of paleobotanical knowledge as compiled in a 
single reference. Accordingly, the data for 
fossil ferns are dealt with here solely for the 
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Taxonomic diversification of ferns based on data compiled by Boureau (1970). Stratigraphic 


resolution of species, genus, and family occurrences from the source is coarser than that of species diversity 
plotted in Figure 1. Consequently, the time-scale of this figure is less precise than desirable. 


purpose of comparing the qualitative diversity 
measurements made at different taxonomic 
levels within a major plant group. 

Figure 2 shows the total diversity of fossil 
ferns based on species, genus, and family 
compilations. Since the citations given in the 
data sets are rarely more specific with regard 
to stratigraphy than lower or upper, middle, 
and lower divisions of each period, the res- 
olution of the geologic scale is necessarily 
more crude than desirable and not directly 
comparable to those given in Figure 1. 

Peaks in the species diversity of ferns occur 
in the Lower and Upper Carboniferous, the 
Middle Jurassic, and in the Eocene. The max- 
imum number of fossil fern species is recorded 
in the Eocene. Species numbers significantly 
decrease in the Middle Carboniferous, Lower 
and Middle Triassic, and in the Oligocene. 
Generic tabulations reveal similar peaks and 
valleys in diversity as those seen in the species 
data; however, the overall diversity maximum 
in the Lower Eocene is less pronounced in 
the generic data set than in the species com- 
pilation. Family diversification is muted 
throughout much of the Phanerozoic. None- 
theless, the Triassic decrease in diversity (from 
ten families in the Upper Carboniferous to six 


families in the Middle Triassic) occurs in both 
the species and generic plots. 

What then is the appropriate taxonomic 
level for gauging long-term patterns in di- 
versification? Large-scale changes in tracheo- 
phyte diversification at the species and family 
levels are poorly correlated (Fig. 1). This is 
not surprising given the diversity of lineages 
that comprise the tracheophytes. The diver- 
sification of tracheophyte species may show 
episodes of increase due to the radiation of a 
few families even if the majority of the re- 
maining families dwindle in species numbers 
or go to extinction. Even for a particular 
group of plants, such as the ferns, there ap- 
pears to be only a loose correlation between 
changes over geologic time in species or genus 
and family numbers, presumably for much 
the same reason—*“‘taxonomic stacking.” 

The family may not be the appropriate 
taxonomic level for analysis for many types 
of evolutionary questions, as, for example, 
the identification of mass-extinction events 
(see next section). A family is eliminated from 
the fossil record with the extinction of its last 
surviving species. This terminal extinction 
event may occur well after the major period 
of species-attrition within a family. Thus, there 
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is a double-edged effect in the use of higher 
taxa for determining unusual perturbations in 
overall diversity. A “‘mass-extinction event” 
may be due to the geologically synchronous 
elimination of the relic species of numerous, 
species-depauperate families. On the other 
hand, many species-rich families could 
undergo a geologically sudden and severe de- 
pauperation in species without the total dis- 
appearance of these families from the fossil 
record. A dramatic reduction in taxonomic 
diversity below the family level would reflect 
a significant event in the history of life that 
could go undetected in family tabulations of 
diversity. Similar arguments could be made 
concerning the effectiveness of genus tabu- 
lations in detecting intense episodes of species 
extinction. 

Species compilations comprise the most bi- 
ologically relevant database for paleontology. 
However, a paleospecies is not operationally 
defined or identified at the same level of taxo- 
nomic scrutiny or with the same body of in- 
formation and experimental format as a species 
of living organism. This severely limits the 
application of conventional species concepts 
to fossil material. 

Recognition of a fossil plant species is per- 
haps even more difficult than that of a fossil 
animal species. Plants are phenotypically more 
plastic than many major groups of animals 
(Cahn & Harper, 1976; Harper & Bell, 1979; 
Harper, 1985; White, 1979). The recogni- 
tion of ecotypes or taxonomic varieties of 
living plants is notoriously difficult and re- 
quires detailed transplant experiments, cy- 
tological and electrophoretic examination, and 
field observations of population dynamics 
across environmental gradients (Harper, 
1977, and references therein). Fossil plants 
are most frequently found as disarticulated 
organs which, for many plants, have been 
previously assigned to separate taxa (Knoll & 
Rothwell, 1981; Stewart, 1983). Only when 
found in organic connection can organ genera 
be placed in synonymy. Consequently, the 
recognition of genus as an “organ” can only 
be done on an ad hoc basis. Since whole plant 
reconstructions are still rare in the paleobo- 
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tanical literature, the inflation of species rich- 
ness due to the effects of dealing with frag- 
mented parts is a serious limitation to assessing 
the true taxonomic diversity in a fossil assem- 
blage. 

Perhaps the only reasonable solution to the 
choice of taxonomic level is to deal with species 
and genus tabulations for each family. Cor- 
relations in the patterns of diversification 
among all three taxonomic levels for partic- 
ular groups of plants could then be used to 
estimate long-term trends or geologically sud- 
den changes in taxonomic richness. As yet 
this suggested type of multiple-level analysis 
has not been undertaken for any group of 
organisms. 


Mass EXTINCTIONS: AN ARTIFACT 
OF HIERARCHY? 


Measurements of diversity based on taxa 
above the species level have been used to 
determine large-scale evolutionary patterns. 
Perhaps the best example of this approach 
comes from the use of family-level data on 
marine animals to determine major diversi- 
fication and extinction patterns (Sepkoski, 
1980, 198la, b, 1984; Raup & Sepkoski, 
1982, 1984). Sepkoski (198 la, b, 1984, and 
elsewhere) has argued that family data rep- 
resent a good proxy for estimating diversity 
of species (cf. Sepkoski et al., 1981). In ad- 
dition, comprehensive data with good strati- 
graphic resolution do not exist for marine 
animal species or even genera (e.g., Sepkoski, 
1984: 247-248). However, as will be shown 
in this section, the behavior of families as 
evolutionary “‘units,”’ and in particular the 
pattern of family extinction rates, is depen- 
dent upon the frequency distribution of species 
within families. 

Raup & Sepkoski (1982, 1984) identified 
five statistically significant “‘mass extinction” 
events in the Phanerozoic record of marine 
animal families: late Ordovician (Ashgillian), 
late Devonian (Frasnian), late Permian (Gua- 
dalupian-Dzhulfian), late Triassic (Norian), and 
late Cretaceous (Maastrichtian). Although the 
magnitudes of species-extinctions are hard to 
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simulation of random species extinction. Data from Cronquist (1968, 1981) were used to construct an initial 
histogram (= iteration 0) representing 57 families and 48,275 species. Random species extinction (20% 
reduction of species number every iteration) results in a “mass species-extinction event” at iteration 1. Sequential 
extinction events (which withdraw fewer and fewer absolute numbers of species with each successive iteration) 


result in shifts of the histogram (cf. Figs. 5, 6). 


estimate from the family data, Raup (1979) 
suggested that the late Permian “mass ex- 
tinction” event involved a reduction of as 
many as 96% of marine animal species. Sig- 
nificantly, Raup’s regression of family back- 
ground extinction rates against time yields a 
negative slope (a long-term decline in family 
extinction rates), which Raup & Sepkoski 
(1982) interpreted as reflecting a ‘“‘fitness- 
optimization” in which surviving marine taxa 
become more resistant to extinction through 
geologic time (but see Kitchell & Pena, 1984). 
More recently, Sepkoski (1984) modeled the 
total rate of background extinction (with con- 
stant extinction parameters) and suggested 
that the decline in rates resulted from differ- 
ential expansion and replacement of evolu- 
tionary faunas. Sepkoski’s model involves lag 
phases in the evolutionary behavior of clades 
within each evolutionary fauna. 

A continuous and constant rate of species 


extinction can produce discontinuous and 
abrupt increases in rates of family extinction 
provided blocks of species go to extinction. 
This can be illustrated by modeling family 
extinction rates based on the species fre- 
quency distributions of angiosperms. 

Data on the approximate numbers of species 
per family in the major groups of angiosperms 
are provided by Cronquist (1968, 1981). His- 
tograms of the species frequency distributions 
in the families of monocotyledons and dicot- 
yledons are shown in Figures 3 and 4. In both 
cases, the majority of families have fewer than 
500 species per family; however, both mono- 
cots and dicots have families with species 
numbers in excess of 5,000. (The Compositae 
with 19,000 species is not shown in the data 
set for dicotyledons.) These distributions are 
similar to those constructed by Clayton (1972, 
1974) for the generic frequency distribution 
in angiosperm families, which conforms to a 
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(cf. Figs. 5, 6). 


logarithmic distribution. Consequently, the 
following arguments against the use of family 
compilations for determining mass extinction 
events can be extended to the use of genus 
tabulations. 

A computer simulation technique was de- 


vised to track the family extinction rate as a 
constant percentage of species is reiteratively 
removed from each clade. The species 
undergoing extinction were selected randomly 
among all families by means of the standard 
linear congruent method, which does not bias 
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extinction (20% reduction in number of species every iteration). Numbers of species decline exponentially; 
numbers of families decline in as near linear function until the 35th iteration. 


the probability of species extinction toward 
large or small families. Consequently, all 
species have an equal probability of going to 
extinction. If a constant percentage of species 
goes to extinction, then a “‘mass-species” ex- 
tinction event occurs on the first iterative 
cycle of each simulation, i.e., each simulation 
starts out with the maximum number of 
species; therefore, the first cadre of species 
to go to extinction is the largest in absolute 
number. Successive extinction events result 
in a progressive decline in the absolute num- 
ber of species undergoing extinction. The use 
of a constant percentage extinction rate pro- 
vides an opportunity to look for possible cor- 
relations between mass-species and mass-fam- 
ily extinction events. Simulations in which a 
constant number of species go extinct (ds/ 
dt = k) rather than a constant percentage of 
species yield similar patterns of family ex- 
tinction but require extended reiterative cycles 
of species extinction. Each simulation based 
on a constant-percentage rate of species ex- 


tinction was run for 50 reiterative cycles. At 
this point most families (99%) became extinct. 

The consequences of successive iterative 
species extinctions on the frequency distri- 
butions of monocot and dicot families is shown 
in Figures 3-5. As expected, each histogram 
gradually shifts toward the left (= numerous 
species-poor families). Figure 5 shows that 
species numbers exponentially decrease and 
asymptotically approach the x-axis (which is 
a very crude analogue to time, since there is 
no a priori way to relate an iteration of species 
extinction to an interval of geologic time); 
decrease in the number of families with suc- 
cessive iterative species extinctions is nearly 
linear. Various species-extinction rates rang- 
ing from 10% to 50% yield comparable pat- 
terns. Significantly, however, total family-ex- 
tinction rates (number of families going to 
extinction per iteration) are not constant. 
Comparisons between the dicot and monocot 
simulations are shown in Figure 6. The total 
family-extinction rate, E,, for dicots increases 
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FIGURE 6. 
dicotyledon and monocotyledon families resulting from 
a 20% species-extinction rate. These values were com- 
puted on the basis of data shown in Figure 5. 


sharply to 7.0 (after 10 iterations), decreases 
to 5.0 (after 15 iterations), and then plateaus 
before it sharply declines after 35 iterations. 
A pattern that is qualitatively and quantita- 
tively dissimilar to dicots is observed for the 
monocot simulation, e.g., Ep never exceeds 
1.4 and increases only slightly over many 
iterations. 

Differences between the patterns of family- 
extinction rates for monocots and dicots are 
due solely to differences in their respective 
initial species frequency distribution in fam- 
ilies. To test the sensitivity of patterns of 
family-extinction rates to the frequency dis- 
tribution of species among families, three ar- 
tificial distributions were created: (1) one in 
which the majority of families contains many 
species (skewed to the right); (2) a distribution 
having an equal number of families in each 
category of species-richness (= horizontal dis- 
tribution); and (3) one in which the majority 
of families contains few species (= the fre- 
quency distribution is skewed to the left). In 
each case, the total number of families was 
held constant (N = 111). Owing to the nature 
of the distributive functions, the total numbers 
of species and families in each distribution 
cannot be equal in each of the three cases. 
However, species-extinction rates are based 
on percentages, and family-extinction rates 
are the subject of concern. Therefore, the 
inequality of species numbers in the three 
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FicuRE 7. Patterns of the total extinction rate com- 
puted for an artificially constructed species frequency 
distribution in which many families are species-rich 
and only a few families are species-poor. The total 
number of families in this simulation is 111; the total 
number of species is 360,871. A “mass species-extinc- 
tion event” is apparent by the fifih iteration of 20% 
species extinction. The largest absolute number of species 
to go extinct occurs at the first iteration. However, data 
are graphed showing the summation of extinction events 
for intervals of five iterations. This is necessary because 
family extinction rates per iteration are very low. A 
“mass family-extinction event” is simulated at the 40th 
iteration. 


frequency distributions does not appreciably 
influence comparisons among the three trends 
in family extinctions. 

Figures 7—9 show the resulting patterns of 
total family-extinction rates for the three ar- 
tificially created frequency distributions. A 
greater resolution of these peaks (in which 
family-extinction events are plotted every it- 
eration) is possible. However, the number of 
families eliminated per iteration is very small 
(1-3 families); therefore, the magnitude of 
peaks in family-extinction rates would be small 
(1-3 families per iteration). By calculating 
family-extinction rates every five iterations, 
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Ficure 8. Patterns of the total extinction rate com- 
puted for an artificially constructed species frequency 
distribution in which all categories of species richness 
are occupied by equivalent numbers of families (N = 
111). The total number of species is 176,600. This 
species frequency distribution produces a “‘spiked” pat- 
tern of total family-extinction rate with a 20% species- 
extinction rate. 


the pattern of family extinctions is empha- 
sized. The distribution skewed to the right 
(= most families with many species) shows a 
pattern vaguely reminiscent of a “mass ex- 
tinction,” in which the rate of family extinc- 
tion varies little about a mean value of 0.12 
through 25 iterations and then sharply in- 
creases to 8.0 after 40 iterations before de- 
clining to 1.6 at 50 iterations (Fig. 7). The 
pattern of extinction occurring up to 25 it- 
erations results from the gradual and random 
extinction of a few species-depauperate fam- 
ilies. With each successive iteration, however, 
species-rich families are reduced in numbers 
of species but, due to their large initial sizes, 
fail to go to extinction. After 25 iterations 
the “‘species buffering” experienced by these 
families is eroded and they become sufficiently 
depleted such that they collectively undergo 
a “mass family extinction.” Only 2.7% of the 
total number of families are eliminated by the 
25th iteration despite the extinction of over 
99.5% of the original number of species. If 
this simulation reflected a real case from the 
fossil record, an abrupt increase in the rate 
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puted for a species frequency distribution with many 
species-poor families and a few species-rich families 
(N = 111). A total of 43,070 species are represented 
in this artificially created distribution. 


of family extinction would not coincide with 
a mass species-extinction event: by the 15th 
iteration, over 96% of all species are elimi- 
nated, despite the survival of all families. 
Similar incongruities between the extent to 
which family extinctions reflect or coincide 
with species extinctions are seen in the other 
two simulations. An even distribution in the 
number of families containing the various cat- 
egories of species numbers produces an as- 
cending but spiked pattern of family-extinc- 
tion rates (Fig. 8). A distribution skewed 
toward many species-poor families yields a 
broad, irregularly plateaued pattern of family- 
extinction rates, which is somewhat similar to 
the pattern generated from the dicot data set 
(Fig. 9). These two frequency distributions 
are clearly artificial since neither is encoun- 
tered in real species, genus, or family tabu- 
lations. Nonetheless, they are of interest be- 
cause they indicate that no frequency 
distribution produces a correlation between 
species- and family-extinction patterns. 
Computer simulations do not indicate that 
family-level compilations are necessarily in- 
adequate to gauge “‘mass extinction” events. 
Clearly, there is no reason to assume that 
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species extinctions are randomly distributed 
among all families. However, (1) “spikes” in 
the regression of family-extinction rates 
against geologic time can reflect the syn- 
chronous demise of many species-poor fam- 
ilies belonging to a once robust taxonomic (or 
grade-level) cohort; (2) patterns in family- 
extinction rates are dependent upon the fre- 
quency distribution of species among the fam- 
ilies considered; (3) mass species-extinction 
and mass family-extinction events may not 
coincide; indeed, the vast majority of species 
within a clade may disappear well before the 
event is noticeable at the family level (as in 
the case of Ginkgo); and (4) analyses of mass 
extinctions require information on the chang- 
ing patterns of species distributions within 
suprageneric taxa. Similar arguments can be 
made for the use of generic tabulations to 
detect mass-extinction events (cf. Raup & 
Sepkoski, 1986; see also Hoffman, 1986). 
Similarly, diversity measured solely on the 
basis of species numbers ignores the relative 
abundance of individuals within taxa—an es- 
sential feature to understanding the ecological 
significance of an extinction event. The “‘mass 
extinction” of many species that contribute 
only a fraction of the biomass to a biota has 
a decidedly different effect than one involving 
species contributing many individuals to a 
biota. 


NONTAXONOMIC MEASUREMENTS 
OF DIVERSITY 


A treatment of diversity strictly from the 
perspective of changes in species numbers 
ignores many significant aspects of plant evo- 
lution. Qualitative changes in reproductive 
and vegetative morphology and quantitative 
changes in the numbers of individuals within 
a taxon are not reflected in patterns of chang- 
ing species numbers. For example, the num- 
ber of gymnosperm species during the Juras- 
sic and Cretaceous remained relatively 
constant despite considerable taxonomic turn- 
over within various lineages and significant 
alterations in morphology and anatomy of 
representative species. Ginkgo biloba, as a 
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species, occurs from the Mesozoic to the pres- 
ent, yet the number of individuals significantly 
declined during the Cenozoic. Clearly, eval- 
uations of diversity based on something other 
than species numbers are important in eval- 
uating adaptive evolution and changes in com- 
munity structure (Spicer & Hill, 1979; Phil- 
lips & DiMichele, 1981; DiMichele et al., 
1985; Fisher, 1985). 

Another reason for evaluating nontaxo- 
nomic measurements of diversity is that in- 
ferences on genomic (hence speciation) rates 
of evolution based on rates of morphological 
evolution may be faulty. Schopf et al. (1975) 
argued that differences in the number of 
species among lineages can be the conse- 
quence of dealing with taxa that differ in their 
degree of morphological complexity. If fossil 
species truly represent genomically distinct 
entities (sensu the biological species concept), 
then rates of morphological and genomic evo- 
lution would be highly correlated. These au- 
thors, however, assert that this direct cor- 
respondence can never be proven. Therefore, 
it is advisable to look at rates of morphological 
evolution and rates of paleospecies origination 
separately. 

Fortunately, changes in species numbers 
and within-assemblage species composition, 
and morphological patterns of long-term evo- 
lution can be dealt with separately to visualize 
tracheophyte diversification. As in most cases 
with the use of paleontological data, however, 
quantitative analyses can only be used to con- 
struct qualitative comparisons or generaliza- 
tions. Among these guarded generalizations 
are: 


(1) The taxonomic composition of fossil plant 
assemblages is rarely, if ever, stable. It 
is altered most significantly by the sub- 
stitution of taxa within related lineages 
during periods of relatively constant over- 
all species numbers. These “‘intrataxo- 
nomic” alterations occur over extended 
geological time-scales and are generally 
discernable among lineages that share a 
common mode of sexual reproduction 


(Knoll, 1986). 
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TABLE 1. 
plants in fossil assemblages for representative intervals. 


Mean species numbers, X, of vascular 


Niklas et al., 
1980 Knoll, 1986 
X (number X (number 
of floras) of floras) 
Early Devonian 5 (7) 4.0 (8) 
Late Devonian 8.2 (9) 10.7 (7) 
Early Mississippian 10 (7) 11.9 (15) 
Late Mississippian 22.5 (2) 23.3 (8) 
Late Jurassic 27.5 (15) 30.8 (6) 
Early Cretaceous 21.5 (15) 30.2 (22) 
Late Cretaceous 43 (7) 54.4 (17) 


(2) Global species numbers generally in- 
crease most significantly with the radia- 
tion of plants sharing a novel mode of 
sexual reproduction (Niklas et al., 1980, 
1985). 

At least within the temporal resolution 
permitted by most paleontological studies, 
changes in taxonomic composition and 
large-scale species numbers are tightly 
linked to changes in vegetative morphol- 
ogy and the appearance of diverse growth 
habits (Niklas, 1987). 

The appearance of new reproductive 
modes and diverse morphologies within a 
clade or cohort of clades may not nec- 
essarily result in the ecological displace- 
ment of previously existing species. 
Throughout much of the Paleozoic and 
Mesozoic, taxonomic radiations are as- 
sociated with the exploitation of environ- 
ments not previously occupied (Tiffney, 
1981; Knoll, 1986). 

The competitive advantages conferred 
upon a taxon by novel reproductive or 
vegetative capabilities are most pro- 
nounced early during its taxonomic ra- 
diation (Knoll et al., 1984; Knoll & Nik- 
las, 1987). 


(3) 


(5) 


Each of these generalities is not without 
exception or debatable inference, since each 
is based on a limited number of studies. It is 
more instructive to review the nature of the 
data upon which these statements are based 
and to examine their deficiencies. 
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To date, only two studies have examined 
changes in the mean species numbers of plants 
throughout the Phanerozoic (Niklas et al., 
1980; Knoll, 1986). In both cases, fossil as- 
semblages (= ‘“‘floras’’) were selected from 
similar depositional environments so as to 
minimize the differential effects of sedimen- 
tological factors on preservation. The floras 
were selected from what were inferred to be 
warm climates. Comparison between the mean 
species numbers per flora for representative 
time periods is given in Table 1. Despite dif- 
ferences in the sources of data, both studies 
reported remarkably similar trends. Mean 
species numbers within floras have increased 
significantly at least twice during the last 420 
million years—doubling between Late De- 
vonian and Late Mississippian floras, and once 
again between Early and Late Cretaceous flo- 
ras. (The lack of an objective baseline com- 
parison for the Late Devonian floras precludes 
evaluation of the data from Early and Late 
Devonian floras.) 

Each of the two increases in mean species 
number coincides with an increase in overall 
species diversity (Niklas et al., 1980, 1985). 
This is to be expected, since the latter incor- 
porates data used to compile within-floras 
species numbers. However, each of the two 
increases in mean species numbers correlates 
with major transitions in the taxonomic com- 
position of floras. This is shown elegantly in 
a study by Knoll (1986), who presented a 
unique analysis of the taxonomic composition 
of floras throughout the Phanerozoic. Figure 
10 (redrawn from Knoll’s study) shows two 
significant resortings of suprageneric groups 
within floras. One occurs with the advent of 
seed plants in the Late Devonian and Missis- 
sippian; another reflects the radiation of flow- 
ering plants in the Cretaceous. If the various 
plant lineages within Knoll’s data set are 
grouped according to their principal modes 
of reproduction, then the changes in taxo- 
nomic composition are seen more clearly (Fig. 
11). As in the plots of large-scale species 
numbers (Fig. 1), the early Paleozoic floras 
dominated by pteridophytes are replaced by 
gymnosperm-dominated Mesozoic floras, 
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which in turn are replaced by floras dominated 
by angiosperms. 

However, Figure 11 masks the taxonomic 
restructuring in floras that occurs within each 
category of reproductive mode. The expan- 
sion of gymnosperm species in Mesozoic floras 
occurs at the expense of pteridophytes made 
up of lineages that did not fall off at equal 
rates. In general, the numbers of fern species 
are much less affected than those of lycopods 
and sphenopsids. Similarly, gymnosperm lin- 
eages undergo a more or less continuous in- 
trataxonomic restructuring. 

The evolutionary appearance of novel 
modes of reproduction, for example, seeds 
and flowers, is often associated with alter- 
ations in growth habit (Tiffney, 1981; Tiffney 
& Niklas, 1985). The number of plant fam- 
ilies characterized by cryptogamic (pterido- 
phytes) and phanerogamic reproduction 
(gymnosperms and angiosperms), as well as 
the number of families with principally non- 
arborescent and arborescent growth forms, 
are plotted through time in Figure 12. In this 
figure a number of “cross-overs” are seen 
which can be related to changes in taxonomic 
composition. For example, during the Car- 
boniferous and Permian, arborescence rises 
even though cryptogamic reproduction re- 
mains more common. Although seed plant 
species dominate Mesozoic floras, nonarbo- 
rescence is more common than arborescence 
during the Jurassic. During the Cretaceous, 
the number of phanerogamic plant families 
gradually increases, until by the Late Cre- 
taceous, arborescent/phanerogamic-domi- 
nated floras which persist through the Ter- 
tiary are established. The Permo-Triassic is 
a period of considerable restructuring in both 
the principal modes of reproduction and 
growth habits. Before the Triassic, arbores- 
cent cryptogamic families are dominant, while 
during the Triassic nonarborescent phanero- 
gamic families are common. 

From a much more comprehensive anal- 
ysis, Tiffney & Niklas (1985) concluded that 
the history of clonality in land plants can be 
segregated into three stages: (1) Silurian to 
lower Carboniferous during which clonal lin- 
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eages dominated, but in which arborescence 
appeared as a vegetative correlate with het- 
erospory or the seed habit; (2) Permo-Triassic 
to Cretaceous, in which families of arbores- 
cent gymnosperms gradually increased and 
gained numerical dominance over families 
of arborescent and nonarborescent pterido- 
phytes; and (3) Cretaceous to present, which 
marks the combination of rhizomatous growth 
and the seed habit (herbaceous angiosperms) 
which became increasingly more important in 
the later Tertiary and Quaternary. 

Unfortunately these analyses are based on 
familial data (whole-plant reconstructions are 
too few to determine large-scale patterns in 
the evolution of tracheophyte growth habits). 
Accordingly they provide few insights into 
species patterns of growth habit or values of 
relative abundance. Although nonarbores- 
cent families dominate much of the early and 
middle Mesozoic, the abundance of arbores- 
cent gymnosperm species in these floras is 
much higher than that of pteridophytes (Fig. 
10). Clearly the family-diversity plots shown 
in Figure 12 are not reflective of community 
structure. 

As Schopf et al. (1975) pointed out, it is 
possible that paleontologists fail to recognize 
the true diversity of morphologically simple 
organisms and overestimate the diversity of 
morphologically complex organisms. Thus, 
periods of rapid taxonomic diversification may 
be inherently related to (and possibly the 
product of) episodes of rapid morphological 
diversification. New vegetative and reproduc- 
tive features provide the potential to discrim- 
inate new phenetic taxa (= paleospecies). As 
the number of potential taxonomic characters 
increases, the number of possible permuta- 
tions of characters increases exponentially. 
There is little agreement among specialists as 
to which feature(s) (anatomy, morphology, 
reproduction, or even geologic age) contrib- 
ute(s) to identification of a new species. Cer- 
tainly among disparate taxonomic groups of 
plants, species are recognized on often very 
divergent categories of features: more derived 
taxa, such as the angiosperms, have more 
numerous and potentially more complex fea- 
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(data redrawn from Knoll, 1986; courtesy of A. H. Knoll). These data show similarities with those presented 
by Niklas et al. (1985) for overall patterns in diversification within separate lineages of vascular land plants. 
Early Paleozoic floras were dominated by archaic pteridophytic groups (rhyniophytes, zosterophyllophytes, 
trimerophytes) which radiated into variously more derived groups during the Carboniferous. Mesozoic floras 
consist of diverse gymnosperm lineages supplemented by numerous ferns. With the radiation of angiosperms, 
the percent taxonomic representation of nonflowering groups declines abruptly during the Late Cretaceous and 


Tertiary. 


tures with which to identify species (e.g., floral 
structure) than more archaic taxa, such as 
pteridophytes (e.g., stelar anatomy, spore-wall 
characters). 

The Devonian flora provides a convenient 
illustration of the potential relationship be- 
tween taxonomic and morphological rates of 
evolution. There is considerable agreement 
among specialists as to the features that dis- 
tinguish Devonian genera and higher taxa. 
Additionally, there are authoritative treat- 
ments of the stratigraphic occurrences of taxa 
and the first and last occurrences of various 
vegetative and reproductive features upon 
which they are based. For example, Chaloner 
& Sheerin (1979) provided a comprehensive 
stratigraphic treatment of Devonian genera 
as well as the first and last occurrences of 


various reproductive and vegetative features. 
Their data are plotted in Figures 13 and 14. 
(Nonvascular genera, such as Sporogonites, 
Parka, Pachytheca, or those having dubious 
status as tracheophytes, such as Taeniocra- 
da, are excluded from these plots.) Anatom- 
ical (tracheid, stelar, and stomatal type; Chal- 
oner & Sheerin, 1979, figs. 2, 3) features 
are plotted separately from reproductive fea- 
tures, such as the position, shape, and type 
of dehiscence of sporangia (Chaloner & 
Sheerin, 1979, fig. 4). The data indicate that 
the number of Devonian genera increases from 
the Pridolian to the Givetian and then under- 
goes a modest decline in the Frasnian and 
Famennian (Fig. 13). By contrast, the number 
of vegetative and reproductive features upon 
which Devonian taxa are based increases 
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FIGURE 11. Percent taxonomic composition of tracheophytes (shown in Fig. 10) converted to represent major 


categories of sexual reproduction. Free-sporing pteridophytes dominate most of the early Paleozoic floras; 
gymnosperms predominate through most of the Mesozoic; and angiosperms dominate the Late Cretaceous and 


Tertiary. 


throughout this interval. Linear regression 
analysis of the total number of taxonomic 
characters versus the number of genera yields 
r = 0.96 (N = 8), which is significant at the 
1% level. This correlation, however, reveals 
very little with regard to the quantitative de- 
crease in Late Devonian genera despite an 
increase in the potential number of taxonomic 
characters in the Frasnian-Famennian. Clear- 
ly, as the number of genera increases it is 
reasonable to expect an increase in the num- 
ber of taxonomically distinguishable features. 
Analyses of the appearance of new genera 
and of new vegetative/reproductive features 
yield poor correlations (Fig. 14). For example, 
a regression of the number of new genera 
against that of the number of reproductive 
features yields r= 0.283. Regressions of 
either the number of vegetative or the total 
number of features against the number of 
Devonian genera yield lower coefficients of 


correlation. Consequently, the taxonomic rec- 
ognition of new Devonian genera does not 
appear to be correlated necessarily with the 
evolutionary appearance of novel reproduc- 
tive or vegetative features. In addition, a re- 
view of the generic descriptions for Devonian 
vascular plants compiled by Gensel & An- 
drews (1984) suggests that paleobotanical 
treatments of early Paleozoic floras are tax- 
onomically conservative. Therefore, at least 
at the generic level, it does not appear that 
estimates of taxonomic diversification in the 
Devonian are artifacts of rapid morphological 
evolution. Since most Devonian genera have 
few species, this conclusion appears warrant- 
ed at the species level as well (cf. Knoll et 
al., 1984). 

Although the perceived taxonomic diver- 
sification of early vascular land plants appears 
not to be biased by rapid morphological evo- 
lution, other episodes of large-scale increases 
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families through the Phanerozoic (redrawn from Tiffney & Niklas, 1985). The scale of the vertical axis changes 
from increments of five families to increments of 20 between the Lower and Upper Cretaceous. 


in numbers of species have not been rigor- 
ously examined. For philosophical and prac- 
tical reasons, therefore, it is reasonable to 
view conservatively the overall patterns of 
tracheophyte species diversification in terms 
of anatomical, morphological, and reproduc- 
tive evolutionary trends, rather than strictly 
in the context of patterns of species origi- 
nation. 


NULL or BIoLoGIcAL HYPOTHESES? 


Historically, evolutionary theory has come 
almost exclusively from observations made on 
living organisms. Charles Darwin was able to 
use the fossil record as evidence for evolution 
but derived his notion of natural selection 
from insights gained from animal and plant 
breeding, biogeography, and natural history. 
Indeed, he found the fossil record singularly 
intractable in supporting many features of his 


theory (Rudwick, 1976). The Modern Syn- 
thesis incorporated paleontology, but even the 
work of George Gaylord Simpson may be 
viewed as an ad hoc rationalization of patterns 
seen in the fossil record based on neontolog- 
ically derived theory (cf. Gould, 1980). Re- 
cently, however, paleontology has generated 
evolutionary hypotheses based on patterns 
seen in the fossil record. This significant shift 
in the source of evolutionary speculation has 
had many effects, not the least of which is a 
re-evaluation of the biases, artifacts, and lim- 
itations that are inherent to paleontological 
data. Clearly, the fossil record can be used 
to generate evolutionary hypotheses, but only 
provided it reasonably reflects biological phe- 
nomena. For a long time paleontologists have 
recognized the numerous geological factors 
that contrive to filter and distort biological 
processes preserved in the record. Much of 
the recent literature focuses on attempts to 
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and vegetative characters used to identify taxa. Data taken from Chaloner & Sheerin (1979, figs. 1-4). Devonian 
genera are exclusive of nonvascular and dubiously vascular plants (e.g., Parka, Pachytheca, Taeniocrada, 
Prototaxites, Sporogonites, Spongiophyton, Protosalvinia): Pri = Pridolian, Ged = Gedinnian, Sieg = Siegenian, 
Ems = Emsian, Eif = Eifelian, Giv = Givetian, Fras = Frasnian, Fam = Famennian. 


deal with these factors and to reconstruct 
information lost from the fossil record (Nich- 
ols & Pollock, 1983). 

The fossil record is the principal source of 
information on long-term patterns of evolu- 
tion. Retained within it are imperfect records 
of taxonomic diversification, major episodes 
of adaptive radiations, and major extinctions. 
This paper has focused on the quality of these 
patterns for vascular land plants. Quantitative 
analyses of the paleobotanical literature pro- 
vide a basis for reconstructing the broad pat- 
terns of floristic and vegetational change oc- 
curring over the last 400 million years. In 
particular, large-scale patterns in numbers of 
species have been used to reconstruct and 
identify two major floristic changes, one at 
the end of the Paleozoic and another toward 
the end of the Mesozoic; within-assemblage 
taxonomic compositions have been used to 
treat broad patterns in vascular plant ecology; 
and morphological/anatomical data have been 
used to reconstruct trends in organography 
and adaptation. As has been seen, however, 


potentially serious limitations and distortions 
exist in the data. In almost all cases, quan- 
titative analyses can be used comfortably only 
to draw qualitative conclusions. 

The most serious difficulty with the fossil 
record comes from attempts to infer mech- 
anisms from patterns. Ecologists are currently 
debating the use of patterns as data, as the 
recent furor over the use of null models attests 
(Harvey et al., 1983). Is there a null hy- 
pothesis for the pattern of species diversifi- 
cation seen in the fossil record? The answer 
is an equivocal yes. The fossil record of land 
plants exists because of processes of non- 
marine clastic and pyroclastic deposition. The 
vast majority of fossil plants are preserved in 
lowland flood plain or lakeside environments 
or are entombed in volcanic ashfalls and mud- 
flows. Consequently, it is conceivable and even 
probable that much of the patterns seen in 
tracheophyte species diversification can be 
explained in terms of factors that influence 
patterns of nonmarine sedimentation. Prin- 
cipal among these is tectonics, which controls 
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FIGURE 14. 
Data from Chaloner & Sheerin (1979, figs. 2—4). 


rates of uplift, erosion, subsidence of sedi- 
ment-accumulating basins, and volcanism 
(Blatt et al., 1972). 

For example, variation in the type of non- 
marine sediments deposited during a geologic 
period could contribute to the apparent pat- 
tern of species diversification. Are periods of 
high numbers of species also those in which 
geological factors favored the deposition of 
sediments in which fossils are easily re- 
covered? Fossil plants are preserved most 
often and with high morphological fidelity in 
fine-grained carbonaceous detrital sediments 
and volcanic ash deposits. Niklas et al. (1980) 
attempted a limited analysis of data on coal 
resources to estimate coarsely the variation 
in nonmarine carbonaceous sediments through 
the Phanerozoic. We concluded that the Car- 
boniferous, Cretaceous, and Tertiary were 
qualitatively different from other geologic pe- 
riods. These three “coal eras” have large coal 
tonnages per unit of outcrop area. High num- 
bers of species for these three coal-rich eras 
could reflect a combination of extensive sedi- 
mentological “sampling” of terrestrial floras, 
excellent preservation (e.g., coal ball petrifac- 
tions), and the consequence of intensive eco- 
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nomic exploration of coal resources. These 
factors probably account for the unusually 
high numbers of species reported for the Car- 
boniferous and Cretaceous-Tertiary. 
Therefore, it is safe to assume that tec- 
tonics and erosional patterns of deposition 
have contributed to the fluctuation in numbers 
of species in substantial ways. Nonetheless, 
these nonbiological or null hypotheses are 
inadequate to explain the patterns of taxo- 
nomic turnover within fossil assemblages or 
broad patterns in plant organographic/repro- 
ductive evolution. Regardless of the quantity 
of sediment deposited or the preservational 
status of plant parts recorded for a geologic 
period, there exists no necessary and suffi- 
cient correlation between the magnitude and 
direction of physical factors operating in the 
fossil record and patterns in taxonomic turn- 
over and morphological evolution. Patterns 
in plant fossil record can be viewed either 
from the perspective of “‘objects”’ (taxa and 
their origination, persistence, and extinction) 
or “properties” (morphological and repro- 
ductive innovations and elaborations). A 
strictly exclusive treatment of either per- 
spective is unlikely to lead to any insights into 
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the relationship between evolutionary pat- 

terns and mechanisms (Sober, 1985). Null 

hypotheses are an essential component to this 
type of inquiry, but the assessment of alter- 
native hypotheses and multiple-causation in 
paleontology requires both geological and bi- 
ological insights. 

The salient conclusions that emerge from 

a review of the plant fossil record can be 

briefly summarized under five points. (1) The 

concepts of a paleospecies and an extant 
species are significantly different (see Gin- 
gerich, 1985). (2) Consequently, the fossil 
record of species diversification is best viewed 
as a document of trends in morphological/ 
reproductive diversification. Quantitative es- 
timates of “‘diversity’’ can be used to draw 
primarily qualitative conclusions. (3) The 
taxonomic richness recorded for any geologic 
period does not provide information on species 
equitability, which is an essential component 
to considerations of paleoecology and the 
ramifications of phenomena such as extinc- 
tion. (4) Although the properties of taxa above 
the level of species are potentially interesting, 
they do not necessarily reflect those of species, 
which are the primary focus of evolutionary 
mechanisms. Finally, (5) much of the pattern 
of taxonomic diversification could be ex- 
plained by factors operating in a strictly geo- 
logical context; however, identification of 
long-term patterns in biological phenomena 
requires the examination of trends in mor- 
phology, anatomy, and reproductive systems, 
in conjunction with patterns of species diver- 

sification (Raup, 1983). 
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